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Abstract 


A parametric study is presented of the buckling behavior of infinitely long sym- 
metrically laminated anisotropic plates subjected to combined loads . The study 
focuses on the interaction of a subcritical (stable) secondary loading state of constant 
magnitude and a primary destabilizing load that is increased in magnitude until buck- 
ling occurs . The loads considered in this report are uniform axial compression , pure 
in-plane bending , transverse tension and compression , and shear . Results are pre- 
sented that were obtained by using a special purpose nondimensional analysis that is 
well suited for parametric studies of clamped and simply supported plates . In particu- 
lar, results are presented for a [±45] 5 graphite-epoxy laminate that is highly anisotro- 
pic and representative of a laminate used for spacecraft applications . In addition , 
generic buckling-design charts are presented for a wide range of nondimensional 
parameters that are applicable to a broad class of laminate constructions. These 
results show the general behavioral trends of specially orthotropic plates and the 
effects of flexural anisotropy on plates subjected to various combined loading condi- 
tions. An important finding of the present study is that the effects of flexural aniso- 
tropy on the buckling resistance of a plate can be significantly more important for 
plates subjected to combined loads than for plates subjected to single-component 
loads. 


Introduction 

Buckling behavior of laminated plates is a topic of 
fundamental importance in the design of aerospace vehi- 
cle structures. Often the sizing of many subcomponents 
of these vehicles is determined by stability constraints, in 
addition to strength and stiffness constraints. One sub- 
component that is of practical importance in structural 
design is the long rectangular plate. These plates com- 
monly appear as subcomponents of stiffened panels that 
are used for wing structures. In addition, long plates 
appear as subcomponents of semimonocoque shells that 
are used for fuselage and launch vehicle structures. 
Buckling results for infinitely long plates are important 
because they often provide a useful conservative estimate 
of the behavior of finite-length rectangular plates, and 
they provide information that is useful in explaining the 
behavior of these finite-length plates. Moreover, knowl- 
edge of the behavior of infinitely long plates can provide 
insight into the buckling behavior of more complex 
structures such as stiffened panels. 

An important type of long plate that appears as a 
subcomponent of advanced composite structures is the 
symmetrically laminated plate. In the present paper, the 
term, “symmetrically laminated plate,” refers to plates in 
which every lamina above the plate midplane has a corre- 
sponding lamina located at the same distance below the 
plate midplane, with the same thickness, material proper- 
ties, and fiber orientation. Symmetrically laminated 
plates remain flat during the manufacturing process and 
exhibit flat prebuckling deformation states. These char- 
acteristics and the amenability of these plates to struc- 
tural tailoring provide symmetrically laminated plates 


with a significant potential for reducing structural weight 
of aircraft and launch vehicles. Thus, understanding the 
buckling behavior of symmetrically laminated plates is 
an important part of the search for ways to exploit plate 
orthotropy and anisotropy to reduce structural weight. 

In many practical cases, symmetrically laminated 
plates exhibit specially orthotropic behavior. However, 
in some cases, such as [±45] s laminates, these plates 
exhibit anisotropy in the form of material-induced cou- 
pling between pure bending and twisting deformations. 
This coupling is referred to herein as flexural anisotropy, 
and it generally yields buckling modes that are skewed in 
appearance, as depicted in figure 1. The effects of flex- 
ural orthotropy and flexural anisotropy on the buckling 
behavior of long rectangular plates subjected to single 
and combined loading conditions are becoming better 
understood. For example, recent in-depth parametric 
studies that show the effects of anisotropy on the 
buckling behavior of long plates subjected to compres- 
sion, shear, in-plane bending, and various combinations 
of these loads have been presented in references 1 
through 4. The results presented in these references indi- 
cate that the importance of flexural anisotropy on the 
buckling resistance of long plates varies with the magni- 
tude and type of the combined loading condition. How- 
ever, the extent of the influence of the combined loading 
condition on the importance of neglecting flexural 
anisotropy in a buckling analysis is not well understood. 

The objectives of the present paper are to present 
buckling results for specially orthotropic plates subjected 
to combined loads in terms of a useful nondimensional 
design parameter and to identify the effects of flexural 



anisotropy on the buckling behavior of long symmetri- 
cally laminated plates subjected to combined loads in a 
more explicit manner than has been given in references 1 
and 2. This second objective is accomplished by model- 
ing various combined loads as a primary system of desta- 
bilizing loads with a secondary system of subcritical 
loads instead of using the traditional approach in which 
each component of a combined loading is treated as a 
destabilizing load with a fixed relative magnitude. This 
approach permits the effects of flexural orthotropy, flex- 
ural anisotropy, and combined loading characteristics on 
plate-buckling behavior to be obtained and presented in a 
direct manner. The primary destabilizing loads that are 
considered consist of uniform axial compression, shear, 
and pure in-plane bending loads; the secondary subcriti- 
cal loads that are considered consist of transverse tension 
or compression and shear loads. Results are presented for 
plates with the two opposite long edges clamped or sim- 
ply supported. A number of generic buckling curves that 
are applicable to a wide range of laminate constructions 
are also presented using the nondimensional parameters 
described in references 1, 2, and 5. 

Symbols 


A m> B m 

displacement amplitudes (see 
eq. (22)), in. 

b 

plate width, (see fig. 1), in. 

D n ,D n , D22> D 66 

orthotropic plate-bending 
stiffnesses, in-lb 

°16* D 26 

anisotropic plate-bending 
stiffnesses, in-lb 

E 2 , Gj 2 

lamina moduli, psi 

"O 

K 

III 

nondimensional buckling coeffi- 
cient associated with critical value 
of an eccentric in-plane bending 
load (see eq. (21) and fig. 1(a)) 

K ^ n *yl\r 

nondimensional buckling coeffi- 
cient associated with critical 
value of a uniform shear load 
(see eq. (20) and fig. 1(a)) 

l)cr 

nondimensional buckling coeffi- 
cient associated with critical value 
of a uniform axial compression 
load (see eq. (18) and fig. 1(a)) 

K y = ( n yOcr 

nondimensional buckling coeffi- 
cient associated with critical 
value of a uniform transverse 
compression load (see eq. (19) 
and fig. 1(a)) 


^x\y = 5 = 0* 
K s\y = 8 = 0* 
K b\y = 5 = 0 

compression, shear, and in-plane 
bending buckling coefficients 
defined by equations (18), (20), 
and (21), respectively, in which 
anisotropy is neglected in the 
analysis 

^1* ^2* L 4 

nondimensional load factors 
defined by equations (14) 
through (17), respectively 

c 

n xV n yb n xyb n bl 

nondimensional membrane 
stress resultants of system of 
destabilizing loads defined by 
equations (10) through (13), 
respectively 

c 

n x2< n y2> n xy2' n b2 

nondimensional membrane 
stress resultants of system of sub- 
critical loads defined by equations 
(10) through (13), respectively 

N 

number of terms in series 
representation of out-of-plane 
displacement field at buckling 
(see eq. (22)) 

N b 

intensity of eccentric in-plane 
bending load distribution defined 
by equation (5), lb/in. 

N xc 

intensity of constant-valued ten- 
sion or compression load distribu- 
tion defined by equation (5), lb/in. 

N*N y 

longitudinal, transverse, and 
shear membrane stress resultants, 
respectively (see eqs. (5), (7), 
and (8)), lb/in. 

N< xl’ N yl’Nxyl’N b i 

membrane stress resultants of 
system of destabilizing loads 
(see eqs. (6) through (9)), lb/in. 

N< x2 ’ N y2< N xy2> N bl 

membrane stress resultants of 
system of subcritical loads 
(see eqs. (6) through (9)), lb/in. 

P’ Per 

nondimensional loading parame- 
ter (see eqs. (14) through (17)) 
and corresponding value at buck- 
ling (see eqs. (18) through (21)), 
respectively 


out-of-plane displacement 
field at buckling defined by 
equation (22), in. 

*.:y 

plate coordinate system 
(see fig. 1), in. 
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<*«,, P> Y* 8 nondimensional parameters 

defined by equations (1), (2), 
(3), and (4), respectively 


D 


Y = 


16 


3 1/4 

(®> 22 > 


( 3 ) 


e o* e i 


r\ = y/&, ^ = x/X 
X 

Xlb 

Vl2 

^mOl) 


symbols that define distribution of 
in-plane bending load (see fig. 1 
and eq. (5)) 

nondimensional plate coordinates 

half-wavelength of buckling 
mode (see fig. 1), in. 

buckle aspect ratio (see fig. 1) 

lamina major Poisson’s ratio 

basis functions used to represent 
buckling mode (see eq. (22)) 


Analysis Description 

Often in preparing generic design charts for buckling 
of a single flat plate, a special purpose analysis is pre- 
ferred over a general purpose analysis code, such as a 
finite element code, because of the cost and effort usually 
involved in generating a large number of results with a 
general purpose code. The results presented herein were 
obtained using such a special purpose analysis. The anal- 
ysis details are lengthy, and only a brief description of 
the analysis is presented. 

Symmetrically laminated plates can have many dif- 
ferent constructions because of the variety of material 
systems, fiber orientations, and stacking sequences that 
can be used to construct a laminate. A way of coping 
with the vast diversity of laminate constructions is to use 
convenient nondimensional parameters. The buckling 
analysis that is used in the present paper is based on clas- 
sical plate theory and the classical Rayleigh-Ritz method 
and is derived explicitly in terms of the nondimensional 
parameters defined in references 1, 2, and 5. This 
approach is an effective method of conducting generic 
in-depth parametric studies of buckling behavior in terms 
of the minimum number of independent parameters 
needed to characterize fully the buckling behavior and to 
obtain results that indicate overall trends and the sensitiv- 
ity of the results to changes in the parameters. The non- 
dimensional parameters used in the present paper are 
given by 


a 


oo 



n 

22 


\l/4 

/ 


P = 


Dl2 + 2D 66 

(D u D 22 ) W2 


( 1 ) 

( 2 ) 


3 1/4 w 

(D n Dh) 

where b is the plate width and X is the half-wavelength of 
the buckle pattern of an infinitely long plate (fig. 1). The 
subscripted D-terms are the bending stiffnesses of classi- 
cal laminated plate theory. The parameters and P 
characterize the flexural orthotropy, and the parameters Y 
and 5 characterize the flexural anisotropy. 

The loading combinations included in the analysis 
are uniform biaxial tension and compression, uniform 
shear, and eccentric in-plane bending, as depicted in fig- 
ure 1. The longitudinal stress resultant N x is partitioned 
in the analysis into a uniform tension or compression part 
and a linearly varying part that correspond to eccentric 
in-plane bending loads. This partitioning is given by 

N X = N xc- N b^ o + fei-eoM ( 5 ) 


where N xc denotes the intensity of the constant-valued 
tension or compression part of the load, and the term 
containing N b defines the intensity of the eccentric in- 
plane bending load distribution. The symbols Eq and 
define the distribution of the in-plane bending load, and 
the symbol r\ is the nondimensional coordinate given by 
tj -y/b (fig. 1). 

The analysis is based on a general formulation that 
includes combined destabilizing loads that are propor- 
tional to a positive-valued loading parameter p that is 
increased monotonically until buckling occurs and inde- 
pendent subcritical combined loads that remain fixed at a 
specified load level below the value of the buckling load. 
Herein, use of the term, “subcritical load,” includes 
loads, such as a transverse tension load, that do not cause 
buckling to occur. In practice the subcritical loads are 
applied to a plate first with an intensity that is below the 
intensity that will cause the plate to buckle. With the sub- 
critical loads fixed, the destabilizing loads are applied by 
increasing the magnitude of the loading parameter until 
buckling occurs. This approach permits combined-load 
interaction to be investigated in a direct and convenient 
manner. 

The distinction between the destabilizing and sub- 
critical loading systems is implemented in the buckling 
analysis by partitioning the prebuckling stress resultants 
as follows: 


Nxc = - NC x1 +N< x2 


( 6 ) 
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N y - -N yl +N y 2 

(7) 

N xy — N X y l + N X y 2 

(8) 


(9) 


where the stress resultants with the subscript 1 constitute 
the destabilizing loads, and those with the subscript 2 
constitute the subcritical loads. The sign convention that 
is used for positive values of these stress resultants is 
shown in figure 1 . The normal stress resultants of the 
system of destabilizing loads, N° xl and N yV are defined 
to be positive valued for compression loads. This con- 
vention results in positive eigenvalues being used to indi- 
cate instability caused by compression loads. 

The buckling analysis includes several nondimen- 
sional stress resultants that are associated with equa- 
tions (6) through ( 9 ). These dimensionless stress result- 
ants are given by 


c K/ 

n ■ = 

7 t 2 (Di,D 22 ) 1/2 

(10) 

Ny/ 
n . = — — — 
yj 2 _ 

n &22 

(11) 

n - N ™ b2 

xyj 2 3 1/4 

« (^11^22) 

(12) 

„ _ "<•/ 
bj n 2 {D n D 22 ) 1/2 

( 13 ) 

where the subscript j takes on the values of 1 and 2 . In 
addition, the destabilizing loads are expressed in terms of 
the loading parameter p in the analysis by 

"x 1 = L 1 ~P 

( 14 ) 

n y i = h P 

( 15 ) 

n xy\ = hP 

( 16 ) 

n bl = L 4 P 

( 17 ) 


where L 1 through L 4 are load factors that determine the 
specific form (relative magnitude of the load compo- 
nents) of a given system of destabilizing loads. Typi- 
cally, the dominant load factor is assigned a value of 1 
and all other factors are given as positive or negative 
fractions. 


Nondimensional buckling coefficients are given by 
the values of the dimensionless stress resultants of the 
system of destabilizing loads at the onset of buckling; 
that is. 




K y ~ ^ n y i >cr 

K s = ( n xy \) cr = 




1/2 = L lPcr 




n 2 D. 


~ L 2 Pc 


22 


W - 

2 3 1/4 ^3 Per 

n (£>11^22) 


( 18 ) 

(19) 

( 20 ) 


K b = ( n bO cr 


F - 

2 /rk _ ,1/2 “ L * P * 

n (^11^22) 


( 21 ) 


where p cr is the magnitude of the loading parameter at 
buckling. Positive values of the coefficients K x and K y 
correspond to uniform compression loads, and the coeffi- 
cient K s corresponds to uniform positive shear. The 
direction of a positive shear-stress resultant acting on a 
plate is shown in figure 1 . The coefficient K b corre- 
sponds to the specific in-plane bending-load distribution 
that is defined by the selected values of the parameters Sq 
and Ej. 

The mathematical expression that is used in the vari- 
ational analysis to represent the general off-centered and 
skewed buckle pattern is given by 


N 

W N&' 1) = X (j4 m sin K $ + B m COS ( 22 ) 

m = 1 

where t, = xfk and q = y/b are nondimensional coordi- 
nates, w N is the out-of-plane displacement field, and A m 
and B m are the unknown displacement amplitudes. In 
accordance with the Rayleigh-Ritz method, the basis 
functions O m (r|) are required to satisfy the kinematic 
boundary conditions on the plate edges at T| = 0 and 
T) = 1 . For the simply supported plates the basis functions 
used in the analysis are given by 


® m 0l) = sin [mini] (23) 


for values of m = 1 , 2 , 3 , ..., N. Similarly, for the clamped 
plates, the basis functions are given by 

O m (T|) = cos[(m-l)rtq]-cos[(m+ l)rcr|] (24) 


Algebraic equations that govern the buckling behav- 
ior of infinitely long plates are obtained by substituting 


4 



the series expansion for the buckling mode that is given 
by equation (22) into the second variation of the total 
potential energy and then by computing the integrals 
appearing in die second variation in closed form. The 
resulting equations constitute a generalized eigenvalue 
problem that depends on the aspect ratio of the buckle 
pattern Mb (fig. 1) and the nondimensional parameters 
and nondimensional stress resultants that are defined 
herein. The smallest eigenvalue of the problem corre- 
sponds to buckling and is found by specifying a value 
of Mb and solving the corresponding generalized eigen- 
value problem for its smallest eigenvalue. This process is 
repeated for successive values of Mb until the overall 
smallest eigenvalue is found. 

Results that are obtained by using the analysis 
described herein have been compared with other results 
for isotropic, orthotropic, and anisotropic plates that have 
been obtained by other analysis methods. These compari- 
sons are discussed in references 1 and 2; in every case the 
analysis described herein was found to be in good agree- 
ment with the results obtained from other analyses. 

Results and Discussion 

Results are presented for clamped and simply sup- 
ported plates that are loaded by various combinations of 
axial compression, transverse tension or compression, 
pure in-plane bending, and shear. For loading cases that 
involve shear, a distinction is made between positive and 
negative shear loads whenever flexural anisotropy is 
present. A positive shear load corresponds to the shear 
loads shown in figure 1. No distinction between positive 
and negative pure in-plane bending loads is necessary for 
flexurally anisotropic plates. 

Results are presented first for [±45] 5 flexurally 
anisotropic plates. The importance of flexural anisotropy 
on the buckling behavior of a [±45] s plate subjected to 
uniform uniaxial compression, pure shear, or pure in- 
plane bending loads is well-known and is documented in 
references 1 and 2. Additional results for this laminate 
are presented herein that show that the importance of 
flexural anisotropy on the buckling behavior depends on 
the type of combined loading. This thin laminate is repre- 
sentative of spacecraft structural components and is 
made of a typical graphite-epoxy material with a longitu- 
dinal modulus £j = 127.8 GPa (18.5 x 10 6 psi), a trans- 
verse modulus E 2 = 1 1.0 GPa (1.6 x 10 6 psi), an in-plane 
shear modulus G 12 = 5.7 GPa (0.832 x 10 6 psi), a major 
Poisson’s ratio v 12 = 0.35, and a nominal ply thickness of 
0.127mm (0.005 in.). The fiber orientation angle of a 
lamina is measured from the plate *-axis to the lamina 
major principal axis denoted by the subscript 1. 

Generic results are presented next, in terms of the 
nondimensional parameters described herein, for a range 


of parameters that is applicable to a broad class of lami- 
nate constructions. The range of each nondimensional 
parameter used herein is given by 0.1 < P < 3.0, 
0 < A, < 0.6, and 0 < 8 < 0.6. The results presented in ref- 
erences 1 and 2 indicate that a value of 0.6 for y and 5 
corresponds to a highly anisotropic plate. For isotropic 
plates, P = 1 and y = 8 = 0. Moreover, for plates without 
flexural anisotropy, y = 8 = 0. Values of these parameters 
that correspond to several practical laminates are given in 
references 1 and 2. 

To simplify the presentation of the fundamental 
generic behavioral trends, results are presented only for 
plates in which y and 8 have equal values (e.g., [±45] s 
laminates). However, these behavioral trends are applica- 
ble to laminates that have nearly equal values of y and 8 
such as a [±35/-15] 5 laminate made of the typical 
graphite-epoxy material described herein. For this 
laminate P = 1.95, y = 0.52, and 8 = 0.51. Further- 
more, results showing the effects of a^, or equivalently 
(Z) n /Z) 22 ) 1/4 , on the buckling coefficients are not pre- 
sented since it has been shown in references 1 and 2 that 
variations in this parameter only affect the critical value 
of the buckle aspect ratio Mb and not the buckling coeffi- 
cient. (The buckling coefficient remains constant.) A 
value of (D n /D 22 ) m = 1 was used in all the calculations 
presented herein. For clarity the compression, shear, and 
in-plane-bending buckling coefficients, defined by equa- 
tions (18), (20), and (21), respectively, are expressed 
fts Kx ly = 5 = 0 * K s L _ § _ o » K b \y -5 = 0 when 
generic results are described for whicn flexural aniso- 
tropy is neglected in the analysis. 

Importance of Anisotropy on Behavior of [±4S] s 

Plates 

Results are presented in figure 2 for a clamped 
[±45] s plate subjected to a destabilizing uniform uniaxial 
compression load. In addition, results are presented for a 
combined-loading condition that consists of uniform 
axial compression and either a uniform subcritical trans- 
verse tension load or a compression load. Similar results 
are presented in figure 3 for a corresponding plate in 
which the destabilizing load is a uniform shear load. 
In these figures the minimum value of the loading 
parameter p found by solving the generalized eigen- 
value problem for a given value of Mb is shown for val- 
ues of Q<Mb<2 . Moreover, the magnitude of the 
subcritical transverse load is indicated in the figures by 
the value of the nondimensional stress resultant n y2 that 
is defined by equation (11). A limiting value of the sub- 
critical load given by n y2 = -4 corresponds to a wide- 
column buckling mode. The dashed lines shown in the 
figures correspond to the solutions of the generalized 
eigenvalue problem that are obtained when flexural 
anisotropy is included in the analysis. In contrast, the 
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solid lines correspond to the solutions that are obtained 
when flexural anisotropy is neglected in the analysis (i.e., 
■Dj 6 = -D26 = 0). The overall minimum value of the load- 
ing parameter for each curve corresponding to n y2 = 0 is 
indicated by an unfilled circle, and the minimum for each 
curve that corresponds to nonzero values of n y2 is indi- 
cated by a filled circle. These minimum values of the 
loading parameter correspond to the value of the buck- 
ling coefficient for each curve. The corresponding values 
of Mb are the critical values of the buckle aspect ratio. 

The results presented in figure 2 indicate that the 
buckling coefficient (minimum loading parameter) of a 
plate that is loaded by pure compression only (n y2 = 0) is 
overestimated by approximately 30 percent of the aniso- 
tropic buckling load and that the critical value of the 
buckle aspect ratio Mb is slightly overestimated if the 
flexural anisotropy is neglected in the buckling analysis. 
For a subcritical transverse tension load given by 
n y 2 = 3.5 (87.5 percent of the magnitude for a wide- 
column buckling mode), the results show that the effect 
of the flexural anisotropy becomes slightly less impor- 
tant; that is, the buckling coefficient is overestimated by 
22 percent for this case. Neglecting the anisotropy for 
this case also yields results that slightly overestimate the 
critical value of the buckle aspect ratio. However, for 
a subcritical transverse compression load given by 
n y 2 = -3.5, the buckling coefficient is overestimated by 
76 percent when the flexural anisotropy is neglected in 
the analysis, and the critical value of die buckle aspect 
ratio is slightly underestimated. 

The results shown in figure 3 for a shear-loaded plate 
with a subcritical transverse tension or compression load 
and the results shown in figure 2 indicate that the effect 
of neglecting the flexural anisotropy in the buckling anal- 
ysis is much more pronounced for the shear-loaded plate 
than for the corresponding compression-loaded plate. In 
particular, the results presented in figure 3 indicate that 
results that are obtained by neglecting the flexural aniso- 
tropy in the buckling analysis of a plate that is loaded by 
shear only (n y2 = 0) overestimate the buckling coefficient 
by approximately 97 percent of the anisotropic buckling 
load compared with 30 percent for the corresponding 
compression-loaded plate. For a subcritical transverse 
tension load given by n y2 = 3.5, the results for the shear- 
loaded plate predict that the buckling coefficient is over- 
estimated by approximately 75 percent if flexural aniso- 
tropy is neglected. More significantly, for a subcritical 
transverse compression load given by n y 2 = -3.5, the 
results predict that the shear buckling coefficient is over- 
estimated by approximately 209 percent when the flex- 
ural anisotropy is neglected in the analysis. For each 
subcritical loading case, the critical values of the buckle 
aspect ratio are overestimated when the anisotropy is 
neglected in the buckling analysis. 


Results that are similar to the results presented in 
figures 2 and 3 are presented in figures 4 and 5 for a 
clamped [±45] s plate that is loaded by a subcritical shear 
load and either a destabilizing uniform axial compression 
load or a pure in-plane bending load, respectively. Two 
groups of curves are shown in the figures that correspond 
to plates without a subcritical load and with either a posi- 
tive or negative subcritical shear load with a magnitude 
equal to 75 percent of the corresponding shear buckling 
coefficient K s (eq. (20)). For the cases in which the 
anisotropy is included in the analysis, the buckling coef- 
ficients for positive and negative shear loads are given by 
K s = 6.12 and K s = —17.16, respectively. For the cases in 
which the anisotropy is neglected, the buckling coeffi- 
cients for positive and negative shear loads are given by 
K s = 12.03 and K s = -12.03, respectively. The results 
shown in figures 4 and 5 therefore include the effects of 
neglecting the anisotropy in the calculation of the sub- 
critical load and the actual calculation of the buck- 
ling coefficient. 

The results shown in figures 4 and 5 indicate that 
neglecting the flexural anisotropy in a buckling analysis 
of the plate that is subjected only to uniform axial com- 
pression or pure in-plane bending (n^ = 0) yields a 
solution that overestimates the buckling coefficient by 
approximately 30 percent of the anisotropic buckling 
load, and the critical value of the buckle aspect ratio is 
slightly overestimated. For the subcritical positive shear 
load with = 0.7 5 K s , the results predict that the buck- 

ling coefficient is overestimated by approximately 
84 percent and 66 percent of the anisotropic buckling 
coefficient for the uniform axial compression loads and 
the pure in-plane bending loads, respectively. Moreover, 
for both loading conditions, the critical value of the 
buckle aspect ratio is also slightly overestimated when 
the anisotropy is neglected in the analysis. In contrast, for 
the subcritical negative shear load with = 0.75 K s , 

the results predict that the buckling coefficient is under- 
estimated by approximately 26 percent and 16 percent of 
the anisotropic buckling coefficient for the uniform axial 
compression and pure in-plane bending loading condi- 
tions, respectively. In addition, the critical value of the 
buckle aspect ratio is slightly underestimated for the uni- 
form axial compression load and slightly overestimated 
for the pure in-plane bending load when the flexural 
anisotropy is neglected in the analysis. 

Generic Results for Specially Orthotropic Plates 

A broad range of buckling results for specially ortho- 
tropic plates that are subjected to combined loads and 
that have either clamped or simply supported edges are 
presented in this section. These results are included to 
demonstrate the effects of flexural orthotropy on plate- 
buckling behavior in a general manner and to provide 


6 



design curves that can be used in conjunction with the 
results presented in the next section to obtain buckling 
coefficients of flexurally anisotropic plates. 

Results are presented in figures 6 through 10 that 
show the generic effects of plate flexural orthotropy on 
the buckling coefficients of clamped and simply sup- 
ported plates with y = 8 = 0. The results in figures 6 
through 8 are for plates that are subjected to a subcritical 
transverse tension or compression load and a destabiliz- 
ing uniform axial compression load, pure in-plane bend- 
ing load, and shear load, respectively. Similarly, the 
generic effects of flexural orthotropy on plates that are 
subjected to a subcritical shear load and a uniform axial 
compression load or to a pure in-plane bending load are 
shown in figures 9 and 10, respectively. The results in 
figures 9 and 10 are also applicable to plates that are 
loaded in negative shear since the shear buckling coeffi- 
cients for specially orthotropic plates that are loaded in 
positive or negative shear have the same magnitude. Two 
sets of curves are shown in figures 6 through 10 for val- 
ues of the orthotropic parameter P = 0.5, 1, 1.5, 2, 2.5, 
and 3. The solid curves correspond to results for clamped 
plates and the dashed curves are for simply supported 
plates. These curves show the buckling coefficient as a 
function of the nondimensional transverse load n y2 in fig- 
ures 6 through 8 and as a function of the nondimensional 
shear load n ^ 2 in figures 9 and 10. 

The results presented in figures 6 through 10 indicate 
that the buckling coefficients of specially orthotropic 
plates (y = 5 = 0) increase substantially as the orthotropic 
parameter P increases. Furthermore, the results in fig- 
ures 6 through 8 indicate that as the subcritical transverse 
load increases through positive values (increasing ten- 
sion), the buckling coefficients increase substantially. 
This trend is shown to be more pronounced for the shear- 
loaded plates than for the compression-loaded plates or 
for the plates that are loaded by pure in-plane bending. 
Moreover, the increase in buckling coefficient with 
increasing subcritical tension load is predicted to be 
slightly more pronounced for the simply supported plates 
than for the clamped plates. The results presented in fig- 
ures 9 and 10 for the plates with subcritical shear loads 
indicate that the buckling coefficient is substantially 
reduced as the magnitude of the nondimensional shear 
load increases. In addition, the results also predict this 
trend to be slightly more pronounced for the simply sup- 
ported plates than for the clamped plates. 

Interaction of Flexural Anisotropy and Loading 

Results are presented in this section that show the 
generic effects of plate flexural anisotropy and combined 
loading condition on the buckling coefficients of 
clamped and simply supported plates. First, results are 


presented for plates that are subjected to a subcritical 
transverse tension or compression load and destabilizing 
uniform axial compression, pure in-plane bending, and 
shear loads, respectively. Then, results are presented for 
plates subjected to a subcritical shear load and to destabi- 
lizing uniform uniaxial compression and pure in-plane 
bending loads, respectively. 

Plates subjected to subcritical transverse tension or 
compression loads . Results are presented in figures 1 1 
and 12, 13 and 14, and 15 and 16, respectively, for plates 
subjected to a subcritical transverse tension or compres- 
sion load and destabilizing uniform axial compression, 
pure in-plane bending, and shear loads. In figures 11 
through 16, the ratio of the anisotropic buckling coeffi- 
cient to the corresponding orthotropic buckling coeffi- 
cient that is computed with y = 8 = 0 (figs. 6 through 10) 
is given as a function of the orthotropic parameter P for 
equal values of the anisotropic parameters (y = 8) rang- 
ing from 0.1 to 0.6. For each value of y = 8 that is given 
in figures 11 through 16, three curves are presented. The 
solid lines correspond to a value of the nondimensional 
transverse load n y2 = 0 that is used in the calculations (no 
subcritical load). Similarly, the finely dashed and 
coarsely dashed lines correspond to values of n y2 = -0.5 
(compression) and 0.5 (tension), respectively, for the 
simply supported plates, and to values of n y2 = -2 (com- 
pression) and 2 (tension), respectively, for the clamped 
plates. The magnitudes of these values correspond to 
50 percent of K y , the buckling coefficient that corre- 
sponds to a wide-column buckling mode (K y = 1 and 4 
for simply supported and clamped plates, respectively). 

The results that are presented in figures 11 and 12, 
respectively, for simply supported and clamped plates 
subjected to a subcritical transverse tension or compres- 
sion load and to destabilizing uniform uniaxial compres- 
sion show that the anisotropic buckling coefficient is 
always less than the corresponding orthotropic buckling 
coefficient for the full range of parameters considered. In 
addition, these results predict that the effects of neglect- 
ing anisotropy are more pronounced for the plates with 
the transverse compression loads than for the plates with- 
out a subcritical load ( n y2 = 0) or a transverse tension 
load. Moreover, these results predict that this trend is 
slightly more pronounced for the clamped plates than for 
the simply supported plates, as indicated by the increase 
in the relative spacing of the three curves that correspond 
to each fixed value of y and 8. The results also predict 
that the reduction in buckling coefficient caused by 
anisotropy is generally larger for the clamped plates with 
n y 2 = “2 (transverse compression at a magnitude of 
50 percent of the wide-column buckling coefficient) than 
for the corresponding simply supported plates. This trend 
is indicated by the observation that the curve for the 
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clamped plates generally lies below the corresponding 
curve for the simply supported plates for each fixed value 
of y and 8. Based on a similar observation, the results 
also predict, however, that the reduction in buckling 
coefficient caused by anisotropy is generally smaller for 
the clamped plates with n y2 = 0 or 2 (transverse tension 
at a magnitude of 50 percent of the wide-column buck- 
ling coefficient) than for the corresponding simply sup- 
ported plates. 

The results presented in figures 13 and 14 for simply 
supported and clamped plates, respectively, that are sub- 
jected to a subcritical transverse tension or compression 
load and a destabilizing pure in-plane bending load, also 
show that the anisotropic buckling coefficient is always 
less than the corresponding orthotropic buckling coeffi- 
cient for the full range of parameters being considered. 
These results also predict that the effects of neglecting 
anisotropy are more pronounced for plates with a trans- 
verse compression load than for plates with n y2 = 0 or the 
transverse tension load, but not to the extent that is 
exhibited by the corresponding plates with the uniaxial 
compression load. (Compare the difference in the rela- 
tive spacing of the three curves that correspond to each 
fixed value of y and 8 that is shown in figures 11 
through 14.) Moreover, these results predict that this 
trend is slightly more pronounced for clamped plates 
than for simply supported plates, but also not to the 
extent that is predicted for the corresponding plates with 
a uniaxial compression load. Furthermore, like the plates 
with a uniaxial compression load, the results also predict 
that the reduction in buckling coefficient caused by 
anisotropy is generally larger for clamped plates with 
n y 2 = -2 (transverse compression) than for the corre- 
sponding simply supported plates. Similarly, the results 
also predict that the reduction in buckling coefficient 
caused by anisotropy is generally smaller for clamped 
plates with n y2 = 0 or 2 (transverse tension) than it is for 
the corresponding simply supported plates. 

The results for simply supported plates subjected to a 
subcritical transverse tension or compression load and a 
positive destabilizing shear load presented in figure 15 
and corresponding results for clamped plates presented in 
figure 16 also show that the anisotropic buckling coeffi- 
cient is always less than the corresponding orthotropic 
buckling coefficient. However, these results also show 
that this trend is reversed for negative shear loads. These 
results also predict that the effects of neglecting aniso- 
tropy are generally more pronounced for plates with 
transverse compression loads than they are for plates 
with n y2 = 0 or with transverse tension loads. This trend 
for the shear-loaded plates is slightly less pronounced 
than it is for the corresponding plates with the desta- 
bilizing uniaxial compression loads and slightly more 
pronounced for the corresponding plates with the desta- 


bilizing pure in-plane bending loads. Moreover, these 
results predict that this trend is also slightly more pro- 
nounced for clamped plates than for simply supported 
plates. Furthermore, like the plates with the uniaxial 
compression loads, the results predict that the reduction 
in buckling coefficient caused by anisotropy is generally 
larger for clamped plates loaded in positive shear and 
with n y 2 = -2 (transverse compression) than for the cor- 
responding simply supported plates. Similarly, the results 
also predict that the reduction in buckling coefficient 
caused by anisotropy is generally smaller for the clamped 
plates loaded in positive shear and with ti y2 = 0 or 2 
(transverse tension) than for the corresponding simply 
supported plates. For the plates loaded by negative shear, 
these trends are reversed. Overall, the results presented in 
figures 1 1 through 16 indicate that plates loaded by posi- 
tive shear typically exhibit substantially larger reductions 
in the buckling coefficient when flexural anisotropy is 
included in the analysis than do the corresponding plates 
loaded by uniform axial compression or pure in-plane 
bending. In contrast, the results for plates that are loaded 
by negative shear predict sizable increases in the buck- 
ling coefficients when the flexural anisotropy is included 
in the calculations. 

Plates subjected to subcritical shear load. Results 
are presented in figures 17 through 20 and in figures 21 
through 24 for plates that are subjected to a subcritical 
shear load and to destabilizing uniform axial compres- 
sion and pure in-plane bending, respectively. Moreover, 
the results in figures 17, 18, 21, and 22 are for a positive 
subcritical shear load, and the results in figures 19, 20, 
23, and 24 are for a negative subcritical shear load. In 
figures 17 through 24, the ratio of the anisotropic buck- 
ling coefficient to the corresponding orthotropic buckling 
coefficient computed with y = 8 = 0 (figs. 6 through 10) 
is given as a function of the orthotropic parameter p for 
equal values of the anisotropic parameters (y = 8) rang- 
ing from 0.1 to 0.6. For each value of y = 8 that is given 
in figures 17 through 24, four curves are presented. The 
solid lines correspond to the value of the nondimensional 
shear load = 0 (no subcritical load). In addition, the 
finely dashed, moderately dashed, and coarsely dashed 
lines correspond to values of n xy2 = 0.25K s , 0.5 K s , 
and 0.8 K s , respectively, where K s is the shear buckling 
coefficient. More specifically, K s is the positive shear 
buckling coefficient in figures 17, 18, 21, and 22, and the 
negative shear buckling coefficient in figures 19, 20, 23, 
and 24. For each case, the shear buckling coefficient is a 
function of p, y, and 8. 

Points on the curves shown in figures 17 through 24 
are obtained by calculating the shear buckling coefficient 
first, assigning the appropriate subcritical value of n^, 
and then computing the buckling coefficient of the 


8 



system of destabilizing loads. Since the calculation of the 
shear buckling coefficient that is used in specifying the 
subcritical loading magnitude also depends on y and 5, 
the results shown in figures 17 through 24 represent the 
overall effect of neglecting flexural anisotropy in both 
the calculation of the subcritical shear loading magnitude 
and the buckling coefficient of the system of destabiliz- 
ing loads. This computational procedure is described 
subsequently in detail by using the results shown in fig- 
ures 25 and 26 for the case of n = 0.8^ and for plates 
with P = 3 that are subjected to a subcritical shear load 
and to uniform uniaxial compression and pure in-plane 
bending, respectively. 

Buckling interaction curves for plates that are sub- 
jected to a subcritical shear load and uniform uniaxial 
compression and to a subcritical shear load and pure in- 
plane bending are shown in figures 25 and 26, respec- 
tively. The solid curve shown in each figure corresponds 
to results for specially orthotropic plates with y = 8 = 0, 
and the dashed lines correspond to results for y = 8 = 0.6. 
Furthermore, the unfilled and filled circles on the curves 
correspond to results for n = 0.8^ and for positive 
and negative shear loading, respectively. The ordi- 
nates of the filled and unfilled circles on the curves 
for y=5 = 0 correspond to the buckling coefficients 
K X L = § = o and ^|y = 5 = o *at are indicated in fig- 
ures 17 through 20 and in figures 21 through 24, respec- 
tively. Similarly, the ordinates of the filled and unfilled 
circles on the curves for y = 5 = 0.6 correspond to the 
buckling coefficients K x and K b that are also indicated in 
figures 17 through 20 and in fig-ures21 through 24, 
respectively. The buckling coefficient ratios shown in 
figures 17 through 24 for the case n xy2 = ^^K s are 
obtained by dividing the value of the ordinate on the 
curve for y = 8 = 0.6 (in figs. 25 and 26) by the corre- 
sponding value of the ordinate for the curve for y = 8 = 0 
for either the filled or unfilled circles. For positive sub- 
critical shear loading, the results in figures 25 and 26 
indicate that the buckling coefficient ratios for plates 
subjected to uniform uniaxial compression and pure in- 
plane bending have values less than 1. In contrast, for 
negative shear loading, the results in figures 25 and 26 
indicate that the buckling coefficient ratios have values 
greater than 1. This trend reversal is caused by the hori- 
zontal shift in the buckling interaction curves that occurs 
when anisotropy is present. 

The results presented in figures 17 and 18 for simply 
supported and clamped plates, respectively, loaded by 
uniform axial compression and a positive subcritical 
shear load, predict monotonic reductions in buckling 
coefficient with increasing values of the anisotropic 
parameters. This same trend is also predicted for increas- 
ing magnitudes of the positive shear load as indicated by 
the four different line types shown in figures 17 and 18. 


Moreover, these results show that the reductions in buck- 
ling coefficients are generally more pronounced for the 
simply supported plates than for the clamped plates. This 
trend is also indicated by the observation that each curve 
for the clamped plates generally lies above the corre- 
sponding curve for the simply supported plates for a 
fixed value of y and 8. More importantly, these results 
indicate that the reduction in buckling coefficient that is 
associated with including anisotropy in the calculations 
increases substantially as the magnitude of the positive 
subcritical shear load increases. 

The results shown in figures 19 and 20 for simply 
supported and clamped plates, respectively, which are 
loaded by destabilizing uniform axial compression and a 
negative subcritical shear load, predict monotonic reduc- 
tions in buckling coefficient with increasing values of the 
anisotropic parameters for the values of 2 = 0 
and 0.2 5K S . However, for the larger values of n^ 2 shown 
in the figures, the results predict increases in buckling 
coefficient with increasing values of the anisotropic 
parameters, as anticipated from the previous results that 
have been presented. This seemingly unusual trend for 
the negative shear load is a manifestation of the horizon- 
tal shift in the compression-shear buckling interaction 
curves caused by flexural anisotropy that has been 
reported in reference 1 and illustrated in figure 25. In 
addition, these results show this trend of increasing buck- 
ling coefficient with increasing values of the anisotropic 
parameters to be, overall, slightly more pronounced for 
clamped plates than for simply supported plates. The 
trend of a decreasing buckling coefficient with decreas- 
ing values of the anisotropic parameters is predicted, 
overall, to be slightly more pronounced for simply sup- 
ported plates than for clamped plates. Furthermore, these 
results indicate that the importance of anisotropy on the 
buckling resistance changes dramatically with the magni- 
tude of a negative subcritical shear load. 

Results similar to those presented in figures 17 
and 18 are presented in figures 21 and 22 for simply sup- 
ported and clamped plates, respectively, that are loaded 
by pure in-plane bending and a positive subcritical shear 
load. These results also predict monotonic reductions in 
buckling coefficient with increasing values of the aniso- 
tropic parameters and increasing magnitudes of a posi- 
tive subcritical shear load. However, these results show 
that the reductions in buckling coefficient are practically 
the same for the simply supported and clamped plates. 
The results also indicate that the reduction in buckling 
coefficient associated with having anisotropy included 
in the calculations increases substantially as the magni- 
tude of the positive subcritical shear load increases. 
Comparing the results in figures 17, 18, 21, and 22 sug- 
gests that the reductions in buckling coefficient that are 
caused by including flexural anisotropy in the analysis 
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are nearly the same for the plates loaded by uniform axial 
compression and pure in-plane bending with ^2 = 0. 
However, the comparison also suggests that the reduc- 
tions in buckling coefficient that are caused by including 
flexural anisotropy in the analysis are more pronounced 
for plates that are loaded by uniform axial compression 
than for those that are loaded by pure in-plane bending 
when a nonzero subcritical shear loading is present, espe- 
cially for the larger magnitudes of a subcritical shear 
load. This trend is also indicated by the observation that 
each curve for the compression-loaded plates generally 
lies below the corresponding curve for the plates that are 
subjected to pure in-plane bending, for a fixed value of y 
and 8. 

Results similar to those presented in figures 19 
and 20 are shown in figures 23 and 24 for simply sup- 
ported and clamped plates, respectively, that are loaded 
by pure in-plane bending and a negative subcritical shear 
load. These results also predict monotonic reductions 
in buckling coefficient with increasing values of the 
anisotropic parameters for the values of n xy2 = 0 and 
0.25£j. Unlike the results for compression-loaded plates, 
the results in figures 23 and 24 predict monotonic reduc- 
tions in buckling coefficient with increasing values 
of the anisotropic parameters for n xy2 = 0.5K s . For 
n xy2 = 0-8ATj, the results predict increases in buckling 
coefficient with increasing values of the anisotropic 
parameters. This unusual trend for a negative shear load 
is also a manifestation of the horizontal shift in the pure 
in-plane bending-shear buckling interaction curves 
caused by flexural anisotropy that has been reported in 
reference 2 and illustrated in figure 26. Moreover, the 
difference between the trends for plates that are loaded 
by negative subcritical shear and pure in-plane bending 
and by negative subcritical shear and uniform axial com- 
pression is manifested by the pronounced difference in 
the shape of the corresponding buckling interaction 
curves shown in figures 25 and 26. (See also refs. 1 
and 2.) For the curves in figures 23 and 24 that show 
increases in buckling coefficient with increasing values 
of the anisotropic parameters, the results predict that the 
increases in buckling coefficient are generally more pro- 
nounced for the clamped plates than for the simply sup- 
ported plates. In contrast, for the curves that show 
decreases in buckling coefficient with increasing values 
of the anisotropic parameters, the results predict that the 
decreases in buckling coefficient are slightly more pro- 
nounced for the simply supported plates tfian for the 
clamped plates. Comparing the results in figures 19, 20, 
23, and 24 indicates that the curves in the figures, which 
show increases in buckling coefficient with increasing 
values of the anisotropic parameters, also show that the 
increases in buckling coefficient are, overall, much more 
pronounced for the compression-loaded plates than for 


the plates loaded by pure in-plane bending. Moreover, 
the results indicate that the curves for plates with nonzero 
subcritical shear loads, which show decreases in buck- 
ling coefficient with increasing values of the anisotropic 
parameters, also show that the decreases in buckling 
coefficient are, overall, much more pronounced for the 
plates loaded by pure in-plane bending than for the 
compression-loaded plates. 

Concluding Remarks 

A parametric study has been presented of the buck- 
ling behavior of infinitely long symmetrically laminated 
anisotropic plates subjected to combined loads. A special 
purpose nondimensional analysis that is well suited for 
parametric studies of clamped and simply supported 
plates has been described, and its key features have been 
discussed. The results presented herein have focused on 
the interaction of a subcritical (stable) secondary loading 
state and a primary destabilizing loading state. The inter- 
action of uniform axial compression, pure in-plane bend- 
ing, transverse tension and compression, and shear loads 
with plate flexural anisotropy and orthotropy have been 
examined. In particular, results have been presented 
for [±45] s thin graphite-epoxy laminates that are repre- 
sentative of spacecraft structural components. In addi- 
tion, a number of generic buckling results have been 
presented that are applicable to a broad class of laminate 
constructions and show explicitly the effects of flexural 
orthotropy and flexural anisotropy on plate-buckling 
behavior under combined loads. 

The most significant finding of the present study is 
that the importance of flexural anisotropy on the buck- 
ling behavior of a long plate is strongly dependent on the 
type and magnitude of the combined load that is applied. 
Specifically, the results presented show that significant 
errors (on the order of 100 percent) can be made when 
calculating buckling coefficients of plates subjected to 
combined loads if flexural anisotropy is neglected. Over- 
all, the results presented herein show that the buckling 
coefficients increase significantly as the orthotropic 

^12 ^^66 

parameter P = — increases, and that they 

(D 11 ^ 22 ) 

decrease significantly as the anisotropic parameters 
£> lfi 

7 = T74 and 5 = 7174 with equa 1 

(^11^22) (^11^22) 

values (y = 8) increase for all load combinations that are 
considered, except those involving negative shear loads. 
For the negative shear loading condition, the trend is 
generally reversed. The results that are presented herein 
also show, in many cases, that the effects of plate aniso- 
tropy are more pronounced for clamped plates than for 
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simply supported plates that are subjected to combined 
loads. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
September 22, 1995 
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Figure 1. Sign convention for positive- valued stress resultants (Eq > ej > 0 shown above). 
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Figure 5. Buckling results for clamped [±45] s plates subjected to pure in-plane bending and subcritical shear loads. 
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Figure 6. Effects of orthotropic parameter P on buckling interaction curves i 
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Figure 7. Effects of orthotropic parameter p on buckling interaction curves for specially orthotropic plates (y = 8 = 0) subjected to pure in-plane bending and 
subcritical transverse tension or compression loads. 
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Figure 9. Effects of orthotropic parameter P on buckling interaction curves for specially orthotropic pi 
critical shear loads. 
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Figure 10. Effects of orthotropic parameter P on buckling interaction curves for specially orthotropic plates (y = 5 = 0) subjected to pure in-plane bending and 
subcritical shear loads. 
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Figure 11. Effects of orthotropic parameter p and anisotropic parameters y and 5 on buckling coefficients for simply supported plates subjected to axial 
pression (N x ^) and subcritical transverse tension or compression (N y2 ) loads. 
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Figure 12. Effects of orthotropic parameter P and anisotropic parameters y and 8 on buckling coefficients for clamped plates subjected to axial compression 
w (A^j) and subcritical transverse tension or compression (N y2 ) loads. 
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Figure 13. Effects of orthotropic parameter p and anisotropic parameters y and 8 on buckling coefficients for simply supported plates subjected to pure 
plane bending (N bl ) and subcritical transverse tension or compression (A^ 2 ) loads. 
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Figure 14. Effects of orthotropic parameter P and anisotropic parameters yand 8 on buckling coefficients for clamped plates subjected to pure in-plane bend 
ing (N b j) and subcritical transverse tension or compression (N 2 ) loads. 
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Figure 15. Effects of orthotropic parameter P and anisotropic parameters y and 8 on buck 
(N xyl ) and subcritical transverse tension or compression (N 2 ) loads. 
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Figure 16. Effects of orthotropic parameter P and anisotropic parameters y and 8 on buckling coefficients for clamped plates subjected to shear (N^y) and 
subcritical transverse tension or compression (AT, 9 ) loads. 
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Figure 17. Effects of orthotropic parameter P and anisotropic parameters y and 5 on buckling 
pression (A^j) and positive subcritical shear (n 2 ) loads. 
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Figure 18. Effects of orthotropic parameter P and anisotropic parameters yand 8 on buckling coefficients for clamped plates subjected to axial compression 
(N x i) and positive subcritical shear (n xy2 ) loads. 
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Figure 19. Effects of orthotropic parameter P and anisotropic parameters y and 8 on bucl 
pression (AT j) and negative subcritical shear (n 2 ) loads. 
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Figure 20. Effects of orthotropic parameter p and anisotropic parameters y and 8 on buckling coefficients for clamped plates subjected to axial compression 
(A^l) and negative subcritical shear (n 2 ) loads. 
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Figure 21. Effects of orthotropic parameter P and anisotropic parameters y and 8 on buckl 
plane bending (N b{ ) and positive subcritical shear (n 2 ) loads. 
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Figure 22. Effects of orthotropic parameter P and anisotropic parameters y and 5 on buckling coefficients for clamped plates subjected to pure in-plane bend 
ing (N b j) and positive subcritical shear (n 2 ) loads. 
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Figure 23. Effects of orthotropic parameter P and anisotropic parameters y and 8 on buckling coefficients for simply supported plates subjected to pure in 
plane bending (N b y) and negative subcritical shear (n 2 ) loads. 







Orthotropic parameter, p = 


D 12 + 2D 66 
(Du D 22 ) m 


Figure 24. Effects of orthotropic parameter P and anisotropic parameters y and 8 on buckling coefficients for clamped plates subjected to pure in-plane bend 
ing and negative subcritical shear (n 2 ) loads. 








Figure 25. Buckling interaction curves for simply supported plates subjected to axial compression (A^, ) and subcritical shear (n^) loads (P = 3, y = 8 = 0, 
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